Integrability of the Pairing Hamiltonian 



M.C. Cambiaggio\ A.M.F. Rivas^*and M. Saraceno^ 



February 9, 2008 

oo 

^Departamento de Fisica 
Comision Nacional de Energia Atomica ( CNEA ), 

cn 

^ \ Ave. del Lihertador 8250, 1^29 Buenos Aires, Argentina. 

O . 

. Centra Brasileiro de Pesquisas Fisicas 

0^ 



Rua Xavier Sigaud 150, CEP 22290-180, RJ, Rio de Janeiro, Brazil 



PACS: 21.69n 21.60Jz 74.20Fg 05.45+b 03.65-w 46.90+s 



^ ' Keywords: Pairing Interaction, Integrability, Time dependent Hartree Fock, Constants of 



the motion, Poincare section. 

Abstract 

We show that a many-body Hamiltonian that corresponds to a system of fermions inter- 
acting through a pairing force is an integrable problem , i.e. it has as many constants of the 
motion as degrees of freedom. At the classical level this implies that the Time-dependent 
Hartree-Fock- Bogoliubov dynamics is integrable and at the quantum level that there are 
conserved operators of two-body character which reduce to the number operators when the 
pairing strength vanishes. We display these operators explicitly and study in detail the 
three-level example. 
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1 Introduction 



Two kinds of simple models are commonly used in nuclear physics for displaying the essential 
properties of the nuclear interaction, the particle-hole and the particle-particle models. The 
simplest one, the Lipkin model[0, consists of two levels of equal degeneracy and fermions 
interacting through a particle-hole force. It has been extensively used to test various approx- 
imation schemes and its classical version, provided by the Time-Dependent Hartree -Fock 
(TDHF) approach, is integrable. When this model is extended to three or more levels ^ Q 
one finds that the TDHF approximation yields a classical problem which is non integrable, 
displays various degrees of chaotic behavior and has an intricate and interesting phase-space 
structure 

In principle the analogous situation for the particle-particle interaction can be thought 
to behave in a similar way. A two-level model with a pairing force is integrable |^ 
and one would expect that the extension to three or more levels would yield non-integrable 
TDHF-Bogoliubov (TDHFB) dynamics. In this paper we report the fact that this is not so 
and that the problem of a pairing force acting in a restricted shell model space with L single 
-particle levels turns out to be integrable both classically and quantum mechanically. We 
display exphcitly the constants of the motion involved and we study their properties, their 
group structure and their classical limits. 

The outline of the paper is as follows. In Section 2 we review the pairing model, in both 
its quantum and classical aspects. The classical limit is obtained by a large scale degeneracy 
argument that leads to the TDHFB equations of motion. Section 3 is devoted to the search 
for new constants of the motion, i.e. new operators commuting with the Hamiltonian. These 
constants of the motion, which are not unique, turn out to be non trivial two body operators 
involving the coupling constant and the single particle energies. A set of commuting operators 
is thus constructed that renders the problem integrable. In Section 4 we treat the three- 



level case and show explicitly the consequences of this integrability both at the quantum and 
classical levels. The last Section is devoted to conclusions and final remarks. 



2 The model 

The pairing force is a very general interacting mechanism that has an ubiquitous role in 
the quantum many body problem. In electron systems it leads to the superconducting 
mechanism and in nuclear physics to the collectivity associated to the pairing degree of 
freedom. In this latter case it provides a simplified description of the short range part 
of the nuclear interaction p]. A schematic model that incorporates this basic mechanism 
can be defined by interacting fermions that can occupy L different single-particle shells of 
degeneracy 2VLi and single-particle energies e^. The fermions interact via a monopole pairing 
force. The Hamiltonian of such a system is 

H = E2q^° - ^{K^K- + K-K^), (1) 
i=i ^ 



where 



K'.=\Y.i^UHm.-\). (2) 

771 7* 
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i=l 4=1 ^ m, 



and 



K- = (K+)t = f]K-. (4) 

i=l 

The {bj-mii^jim^) are the usual fermion operators obeying anticommutation relations which 
create or annihilate a fermion on the i-th shell which has degeneracy 2ji + 1 = 2Qi. The 
operators count pairs of particles Ni in each shell by Kf = {Ni — . The operators 
{K^ , iTj", Kf) conform an SU{2) algebra whose Casimir operator is 

k^ = Kf + ^{KtKr + KrK+). (5) 



The full dynamics of the system occurs in the group space of [SU{2)]^ x [SU{2)]2 x . . . x 

[sumL- 

The classical limit of the model is obtained from the TDHFB approximation when the 
degeneracy fij of each level goes to infinity 0. ^l^^ is the semiclassical parameter analogous 
to h in the usual semiclassical treatments. 

One way of obtaining this limit is through the time dependent variational principle imple- 
mented through coherent states that are constructed from the vacuum (or minimal weight) 
state |0), characterized by 

Kr|0) = OK°|0) = -^|0) (6) 



The coherent state in this representation is [|I0| 

\z) = \z^...ZL) = e^ti'^^'\0). (7) 

The equations of motion obtained through the time-dependent variational principle with 
this state are equivalent to the TDHFB equations. To obtain them we use the variational 
principle appropriate for non- normalized states B with an action defined as (we set h = 1) 



S= dt 



2 {^\^P) {^P\^) 
A detailed derivation is given in [Q. The variables Zi are not canonical but the transformation 



u!i = \ -^Zi with < UitUi < (9) 

V 1 + ZiZi 

yields canonical variables satisfying {ui,ujj} = 6ij. The finite range of these variables is a 
consequence of the Pauli exclusion principle between correlated fermion pairs. In terms of 
LJi the variational equations become ordinary hamiltonian equations in complex form, 



where Ti = liniQ-^oo ^^j/^^^ is the classical hamiltonian associated to the problem. In the 
same way for any quantum operator A its classical limit in this representation is, 



Therefore the operators {Kf^K^ and Kf) have their classical analogues, (/C^,/Cj , /C! 
which can be written in terms of the Ui variables 

/C° = uj^Ui-^ (12) 



K-i = ujiJVti- ujiUJi (13) 



/Cj = uJi\IVti - UiUJi. (14) 

This last set of operators obey the classical SU{2) Poisson bracket relations [Q. 

The classical pairing hamiltonian can then be expressed, in analogy with the quantum 
one, directly in terms of the SU (2) generators, 

'H = j:2eX^-G{lC^K- + j:%. (15) 

i=i i=i 



Energy conservation is guaranteed by the time-dependent variational principle |Tl| so that 
the motion occurs in the (2L — l)-dimensional manifold defined hy 1-L{uj,uj) = E. There is a 
further constant of the motion linked to the conservation of the total number of pairs ||^ 

M{uj,uj) = J2 ^^l^'l^^ = Y.^i'^i- (16) 

i=i < > i=i 

The TDHFB equations are then classical hamiltonian equations in a phase space of 2L 
variables. The conservation of H and J\f implies that the case L = 2 is classically integrable, 
a fact that was exploited in to compute energy levels and transition matrix elements 
semiclassically. 

3 Search for additional constants of the motion 



The proof of the integrability for L > 2 requires the existence of L — 1 independent, well de- 
fined , global functions (constants of the motion), whose Poisson brackets with each other and 



with the hamiltonian vanish. Following Hietarinta[|T2[ a quantum mechanical hamiltonian of 
L degrees of freedom is defined to be quantum integrable if there are L — 1 independent, well 
defined, global operators (quantum invariants) which commute with each other and with 
the hamiltonian. Then the energy spectrum of a quantum integrable hamiltonian system is 
naturally labeled by L quantum numbers , which are the eigenvalues of the corresponding 



quantum invariants. Likewise the stationary states are simultaneous eigenfunctions of the L 
corresponding operators. 

For the cases L = 1 and L = 2 the integrability is trivial. H and Kq (or A^) provide the 
commuting operators and 7i and /Co (or A/") the corresponding classical conserved quantities. 
But in the case of L > 2 no new obvious quantum invariants are present and we could expect 
the system to be non integrable and therefore display generically regions of chaotic behavior. 
However we now show how to construct non trivial operators, independent of the hamiltonian 
and the total number of pairs and commuting with them, which make the problem integrable. 

For this purpose let us construct the more general two-body operator O which is hermitian 
and conserves the total number of pairs 



R = j:ia,,K^K^ + P^.K^KT) + ^ j,K^ (17) 
with aij = aji and (3ij = (3ji, and require it to commute with the hamiltonian 

[H, R] = 0. (18) 

The following relations among its coefficients 

aij = f3ij (19) 
_ (j ag + 7t — ajj — 7j 

~ 2 (e,-6,) ' ^'"^ 
are sufficient for the commutativity but they do not determine completely the coefficients. 
Thus several solutions can be found, and one should further check that the operator con- 
structed is independent of H and A^. 



For example if we choose 



an = Oand7i = 



(21) 



we then would have 
R 

where 



H G 
T ~ 2" 



i=l 



(22) 



(23) 



is the Casimir of the SU{2) algebra in the ith-level. In this case R is not useful because it 
is not independent of the conserved magnitudes that we already know. 
However another choice 



an = and 7^ = 5i, 



leads to a set of operators 



(24) 



(25) 



with ki.kj = K^K^ + \{K+KJ- + KrK+). 

In the non interacting limit {G — 0) these operators become the natural set of commuting 
operators {K^; i = 1, . . . , L}. A straightforward but tedious calculation shows that 



(26) 



for any value of G and e^. On the other hand, it is easy to see that H and N can be written 
in terms of the Ri as 

L 



and 



i=l 



H = Y, 26,i?, + G{Y. Rif -Gj2k- 



(27) 



(28) 



i=l 



1=1 



i=l 



Therefore we have constructed a set of L commuting operators which also commute with 
the pairing hamiltonian and are number conserving. They extend to the fully interacting case 
the trivial properties of the number operators of the G = system. Consideration of these 
operators then demonstrate the integrability of the quantum problem. The simultaneous 
eigenvalue equations 

R,\^{\i...\l)) = M'^{\i...\l)) (29) 
gives the eigenvalues of the hamiltonian as 

H = j2 2e.A, + G{j: - J - !)■ (30) 

1=1 1=1 ^ i=l 

It should be noticed however that the actual solution is by no means simplified by this knowl- 
edge. The operators Ri are two body operators as complicated in principle as the hamiltonian 
itself and they cannot be used (except by diagonalizing them) to separate the hamiltonian 
into invariant subspaces. however, as we will see, the simple fact of their existence has very 
drastic implications on the structure of both eigenvalues and eigenf unctions. 

In the classical limit the associated operators for the set {Rf, i = 1 . . . L} written in terms 
of the classical operators (/C^, /C~, /C°) are 

7^. = /C°-GE^^^ (31) 

where JCiXj = iC^lC] + ^/C+ZCj + /Cr/C+). 

Taking into account (|l]), (0), (p5|) and (|3TD we can see that the classical operators TZi,Af 
and H have the same structure (in terms of the SU (2) algebra generators ) as their quantum 
analogs, except for additive constants. It is then easy to see that 

{7^,7^,} = o {^^,ni} = o (32) 

and {7^^,7^J} = 0. (33) 

(34) 



Analogously to the quantum case, the mean number of pairs M and the energy 7i are 

L 

M{ni...nL) = Y.^i (35) 

i=l 

and 

7^(7^l . . . 7^L) = E 2e,7^, + Gijz n,f -^E^l (36) 

i=l 1=1 i=\ 

The functions TZi{u},uj) are surfaces on the 2L dimensional phase space. The trajectories lie 
in the intersection of these surfaces, which are L dimensional tori labeled by the constants 
TZi- Chaotic motion therefore cannot occur in this system. 



4 Manifestations of integrability in the three-level case 

In this section we restrict our analysis to the three- level case (L = 3) and equal degeneracy 
Q^i — Vl and show the consequences of the integrable behavior, both in the quantum and in 
the classical solutions. 

As in the previous section we will start with the quantum treatment. Although the 
problem has three degrees of freedom we have already seen that the total number of pairs 
N in the system is a conserved magnitude. We can then use this fact to reduce explicitly 
the dimensionality to two degrees of freedom and therefore we only need to display another 
quantum invariant operator to have an integrable quantum problem. We then choose it as 
O2 which is a linear combination of the Ri operators defined in the previous section, 

02 = j2^lRi- (37) 

i=l 

The quantum integrability shows up clearly when we display a grid of the simultaneous 
eigenvalues of H and O2. This is done in Fig. 1 for the case where the total number of 
pairs is N — Q — lA, which corresponds to a set of 120 levels. We find that the eigenvalues 
lie on a regular grid that includes all eigenvalues. The fact that this grid is not parallel 



and equally spaced reflects the fact that O2 and H are not action variables that quantize at 
integer spaced values. Of course a point transformation to a set Si and 5*2 exists if O2 and 
H are independent. But we do not construct this variables explicitly. It is clear however 
that the grid we obtain is a smooth deformation of a regular one. This proves that the two 
operators (O2 and H) are independent and that their simultaneous eigenvalues form a set of 



good quantum numbers] 13 



For the classical analysis we introduce the non-interacting action-angle coordinates, 

Ui = uJiLUi (38) 
(pi = arg{ui) (39) 

where Ui is the mean number of pairs in the level i and is now a continuous classical variable. 
In these variables the classical hamiltonian is, 

n =2elr^l + 2e3n3-^](el + e3)- G {ni(fi - ni) + ^^2(^^ - ns) + ^3(1] - ns) (40) 



^1^n\n2\J (r^ - nx){VL - 77.2) cos(0i - ^2) 



^-2^nxn^\l {VL - n^iVt - nj) cos(0i - 03 



^-2^n2n-i^ {VL - n^) iVt - 77,3) cos(02 - 0; 

where we have taken £2 = as the energy reference. 

The analysis is best performed by explicit elimination of the conserved quantity M 
J2i'>T'i- We introduce the canonical transformation. 



/i = ^ ^1 = 01 - 03 (41) 

^2 = 6*2 = 02 - 03 

^3 = ^(^1 + ^^2 + ?^3) = 1 6^3 = 03. 

We have scaled the variables so that they are all in the same range and adopted the mid-shell 
value, i.e. Af{uj,u;) = Q. The problem is reduced to a space of two degrees of freedom whose 



effective hamiltonian is, 

n{h,l2,0u02) = 2einh + 2e3n{l-h-l2)-n{e, + es) (42) 

- Gn^ - Ji) + hii - h) + (/i + /2)(i - /i - /2) 

+2^h{l -h- /2) V(l - + h) cos(^i) 

+2^/12(1 - Jl - /2)v/(l - /2)(/l + /2) C0S(^2)} . 

In the above coordinates the classical version of O is 

02(/i,/2,^i,^2) = elnh-eln{h + h) + ^{el-el) (43) 

- GO? |ei/i(i - h) + e3(/i + /2 - ^)(1 - /i - /2) 
+ei/r^7(l-/i)(l-/2)cos(ei - ^2) 

+ (ei + £3)^/^(1 - /i - /2)\/(l - + h) cos(^^i) 

+63^/2(1 - /l - /2) V(l - /2)(/l + /2) C0S(^2)} . 

It can be explicitly verified that {02,'H} = 0. We have also tested in the numerical 
integration of the equations of motion that both Ti and O2 are constants; and use this fact 
to control the numerical accuracy of the solutions. Integrability is also apparent in Poincare 
sections, as Fig. 2 shows. Only separatrices (no chaotic layers) are observed for any value 
of the coupling constant. The salient feature of this figure is the appearance of forbidden 
regions, quite a common occurrence in spin systems Q]. 

5 Conclusions 

Our analysis shows that the simple pairing force is integrable both at the quantum and at 
the classical levels. The reason why this simple feature has escaped attention can be found 
in the fact that the pairing force has been used almost exclusively in the context of quantum 



many-body physics and for energies close to the ground state. Thus there are no significant 
consequences for ground state or RPA modes, which are in any case integrable. However 
the structure of the highly excited states and eigenfunctions should show this consequence. 
For example, the fiuctuation properties of the eigenstates will follow a Poisson rather than 



a GOE statistics [|14|, while eigenfuntions will show the traces of conserved quantities and 
will not be ergodic. 



There are very few models of interacting particles that are integrable |]T5[ and the fact 
that this simple model belongs to this class comes as a surprise. We have reported some 
consequences in the TDHFB dynamics and on the spectrum. Other features should follow 
naturally, for example the statistics of levels should be very different than the particle-hole 
case0. Also the addition of a small perturbation to this hamiltonian should follow the 
KAM[|16| perturbation scheme and not a fluctuating spectrum characteristic of the pertur- 
bation of chaotic systems. From the technical side it is quite useful to have at one's disposal 
a model which takes into account important features of the nucleon-nucleon interaction and 
which remains integrable for arbitrary values of its strength. Thus perturbation expansions 
can have a steadying point which need not be the non interacting fermion system. 
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FIG. 1. Simultaneous eigenvalues of the Hamiltonian (H) and O2 for = 14 and G = 0.1. 
The variables have been scaled dividing them by N. 

FIG. 2. Poincare surface of section for the three-level pairing hamiltonian . The mid-shell 
value J\f — fl has been taken and the section is performed in the Ii — case (see text). 
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